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Abstract 

We present a measurement of the D^-D^ mixing parameter ycp using a flavor- untagged sample 
of — > K^K~^K~ decays. The measurement is based on a 673 fb~^ data sample recorded 
with the Belle detector at the KEKB asymmetric-energy e^e~ collider. Using a method based 
on measuring the mean decay time for different K^K~ invariant mass intervals, we find ycp = 
(+0.11 ± 0.61(stat.) ± 0.52(syst.))%. 

PACS numbers: 13.25.Ft, 12.15.Ff 
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I. INTRODUCTION 



Particle-antiparticle mixing has been observed in the neutral kaon, and meson 
systems, and evidence for mixing has recently been found for neutral D mesons. The mixing 
occurs through weak interactions and gives rise to two distinct mass eigenstates: |-Di,2) = 
p\D^) ± g|-D°), where p and q are complex coefficients satisfying |pp + |gp = 1. The time 
evolution of the flavor eigenstates, D° and Z)°, is governed by the mixing parameters x = 
(mi— m2)/r and y = (Fi — r2)/2r, where mi2 and ri^2 are the masses and widths of the two 
mass eigenstates D12, and T = (ri + r2)/2. In the Standard Model (SM) the contribution of 
the box diagram, successfully describing mixing in the B- and /T-meson systems, is strongly 
suppressed for mesons due both to the smallness of the Vub element of the Cabibbo- 
Kobayashi-Maskawa matrix [1], and to the Glashow-Illiopoulos-Maiani mechanism [2]. The 
largest SM predictions for the parameters x and y, which include the impact of long distance 
dynamics, are of order 1% [3]. Observation of large mixing could indicate the contribution 
of new processes and particles. 

Evidence for D^-D^ mixing has been found in — s> K~^K~ /tt^tc^ [4, 5], K^tt~ 
[6, 7] and —>■ K^h'tt^ [8] decays. Currently the most precise individual measurements of 
mixing parameters are those from the relative lifetime difference between decays to CP 
eigenstates and flavor-specific final states, ycp-, which equals the parameter y in the limit 
where CP is conserved. Thus far, only CP-even final states K^K~ and tt^tt' have been 
used; the resulting world average value [9] for ycp is (+1.13 ± 0.27)%. 

In this paper we present a fiavor-untagged measurement of ycp using the CP-odd com- 
ponent of KgK~^K~ decays [10]. The measurement is performed by comparing mean 
decay times for different regions of the three-body phase space distribution. As this method 
does not use a fit to the decay time distribution, it does not require detailed knowledge 
of the resolution function or the time distribution of backgrounds. The result has similar 
statistical sensitivity to that obtained by fitting the decay time distribution. 



II. METHOD 



The time-dependent decay amplitude of an initially produced or can be ex- 
pressed in terms of the neutral D meson amplitudes {K'^K+K-\D°) = A{so,s+) and 
{KlK+K-\D^) = A{so,s+), where sq 

expressions are [12, 13]: 



M|.+^_ [11] and s± 
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with ei^2(t) = exp{— z(mi_2 — iiri^2/2)t}. In the limit of CP conservation {p/q = 1), Eqs. (1) 
and (2) simplify to: 

{KlK+K-\D%t)) = Aiiso, s+)ei(t) + A(so, s+)e2{t) (3) 

(KlK-'K-lWit)) = Aiiso, s+)ei(t) - ^2(^0, s+)e2{t) , (4) 

where A(so, s+) = [A{so, s+) + A{so, s+)]/2 and A(so, s+) = [A{so, s+) - A{so, s+)]/2. In 
the isobar model the amplitudes A and A are written as the sum of intermediate decay 
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FIG. 1: Projection of time-integrated Dalitz distribution to sq (solid line), and the ai(so) (dotted 
line) and a2('So) (dashed line) contributions for the Dalitz model given in [14]. 

channel amplitudes (subscript r) with the same final state, A{so, s+) = J2r are*'^''v4.r(so, s+) 
and A{sq, s+) = a^e^'^'-^j.lso, ■§+) = J2r (ir^^'^''A.r{so, s_), where CP conservation in decay 
has been assumed in the final step. If r is a CP eigenstate, then Ar{sQ, s_) = ±Ar{so, s+), 
where the sign +(— ) holds for a CP-even(-odd) eigenstate. Hence the amplitude Ai is 
CP-even, and the amplitude A2 is CP-odd. 

Upon squaring Eqs. (3) and (4) we obtain for the time-dependent decay rates of initially 
produced and 

^^^^^1^ oc |^i(so,.+)|V^(^+^) + |^,(so,.+)|V^(^-^) 
+2Pe[^i(so, S-^)Al{so, s+)] cos (^~^ 
+2Im[Ai{so, s+)A*2{so, s+)] sin (^^-^^ (5) 
oc |^i(.o,3+)|V^(^+^) + |^,(.o,3+)|V^(^-^) 
— 2Pe[^i(so, s+)Al{so, s+)] cos '^"^ 
-2Im[Ai{so, s+)A*2{so, s+)]sm ^x^^e"-, (6) 

where r = l/F is the lifetime. It can be shown (see Appendix) that in the projection of 
the Dalitz plot onto sq, the last two terms in Eqs. (5) and (6) vanish. Hence, a projection 
onto So of the time-dependent decay rate for — > K^K^K~ in the limit of CP conservation 
depends only on the mixing parameter y: 

oc ai(so)e-^(^^^) + a,(so)e-^(-^) , (7) 

where 01,2(59) = / |v4i,2(so, s+)\^ds+. 

Figure 1 shows the time-integrated projection of the decay rate (Eq. 7) together with 
the ai(so) and 02(^0) contributions; the plots are obtained using the Dalitz model of 
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Ref. [14] and taking y = 0. The Dahtz model includes five CP-even intermediate states 
(J'S:gao(980)0, /s:g/o(1370), i^g/2(1270), /s:gao(1450)°, K^/o(980)), one CP-odd intermediate 
state {Kg(j){1020)) and three fiavor-specific intermediate states (-ft'~ao(980)"^, i^~ao(1450)^, 
K+aol980)-). 

The two terms in Eq. (7) have a different time dependence as well as a different Sq 
dependence (see Fig. 1). In any given sq interval, TZ, and assuming y <^ 1, the effective D'^ 
lifetime is 



the events in the interval 71 due to the Ai amplitude. In Eq. (8) we introduced the usual 
notation ycp for the mixing parameter y to indicate that we assumed CP conservation 
in deriving Eq. (7). The definition of ycp in Eq. (8) is consistent with that used in the 
— s> K^K~ /tt^tt~ measurement [4]. 

The mixing parameter ycp can be determined from the relative difference in the effective 
lifetimes of the two Sq intervals, one around the 0(1020) peak (interval ON) and the other 
in the sideband (interval OFF). Using Eq. (8) and taking into account the fact that [1 — 
(/on + foFF)\ycp < 1, we obtain 



The sizes of the ON and OFF intervals are chosen to minimize the statistical uncertainty on 
ycp- They are determined using the Dalitz model of K^K^K~ decays from Ref. [14]. 

The optimal intervals are found to be: M^+k- ^ [1.015, 1.025] GeV/c^ for the ON interval, 
and the union of Mk+k- e [2mx±, 1.010] GeV/c^ and Mk+k- e [1.033, 1.100] GeV/c^ for 
the OFF interval. 

III. MEASUREMENT 

This section is organized as follows: in subsection III A we describe how signal decays 
are reconstructed; in subsection III B we describe how the mean decay time of the signal is 
extracted in the presence of background; in subsections III C and III D we describe how the 
background fraction and mean lifetime, respectively, are determined; in subsection III E we 
describe how /on — /off is determined; and in subsection IIIF we give the result for ycp. 

A. Reconstruction of events 

The data were recorded with the Belle detector at the KEKB asymmetric-energy e~^e~ 
collider [15]. The Belle detector consists of a silicon vertex detector (SVD), a 50-layer 
central drift chamber (CDC), an array of aerogel threshold Cherenkov counters (ACC), a 
barrel-like arrangement of time-of-flight scintillation counters (TOF), and an electromagnetic 
calorimeter (ECL) comprised of CsI(Tl) crystals located inside a superconducting solenoid 
coil that provides a 1.5 T magnetic field. An iron fiux-return located outside of the coil is 
instrumented to detect mesons and to identify muons (KLM). The detector is described 
in detail elsewhere [16]. Two inner detector configurations were used. A 2.0 cm beampipe 
and a 3- layer silicon vertex detector was used for the first sample of 156 fb~^, while a 1.5 





(9) 



6 



cm beampipe, a 4-layer silicon detector and a small-cell inner drift chamber were used to 
record the remaining 517 fb~^ of data [17]. We use an EvtGen- [18] and GEANT-based [19] 
Monte Carlo (MC) simulated sample, in which the number of reconstructed events is about 
three times larger than in the data sample, to study the detector response. 

The candidates are reconstructed in the tt^ti~ final state. We require that the pion 
candidates form a common vertex with a fit probability of at least 10~^, and that they be 
displaced from the e"'"e~ interaction point (IP) by at least 0.9 mm in the plane perpendicular 
to the beam axis. We also require that they have an invariant mass M^+^- in the interval 
[0.468, 0.526] GeV/c^. We reconstruct D° candidates by combining the candidate with 
two oppositely charged tracks assumed to be kaons. We require charged kaon candidate 
tracks to satisfy particle identification criteria based upon dE/dx ionization energy loss in 
the CDC, time-of- flight, and Cherenkov light yield in the ACC [20]. These tracks are required 
to have at least one SVD hit in both r — (p and z coordinates. A momentum greater 
than 2.55 GeV/c in the e+e~ CM frame is required to reject D mesons produced in S-meson 
decays and to suppress combinatorial background. Events with a KgK~^K~ invariant mass 
(M^o^+^-) in the interval [1.77, 1.96] GeV/c^ are retained for further analysis.. 

The proper decay time of the D° candidate is calculated by projecting the vector joining 
the production and decay vertices, L, onto the momentum vector pjj: t = {m£,o/pu)L ■ 
[Pd/pd), where m^jo is the nominal mass. Charged and neutral kaon candidates are 
required to originate from a common vertex for which the fit probability is larger than 
10~^. According to simulation studies, if the decay position is determined by fitting the 
two prompt charged tracks to a common vertex, the decay length and the opening angle of 
the and K~ (and thus their invariant mass) are strongly correlated. This correlation 
is avoided by determining the decay length from a fit where only a single charged kaon 
and the K'^ are fitted to a common vertex. Both K^Kg vertex combinations are required 
to have a probability larger than 10~^; for the L determination, the one with the higher 

fit probability is chosen. The D'^ production point is taken to be the intersection of 
the trajectory of the candidate with the IP region. The average position of the IP is 
calculated for every ten thousand events from the primary vertex distribution of hadronic 
events. The size of the IP region is typically 3.5 mm in the direction of the beam, 100 fim 
in the horizontal direction, and 5 /im in the vertical direction. The uncertainty in a D^'s 
candidate's proper decay time (aj is evaluated from the corresponding covariance matrices. 
We require at < 600 fs. The maximum of the at distribution is at ~ 230 fs. 

Around 362 x 10^ events pass all selection criteria. The {Mj^+t^-, M^oj^+^-) and M^+k- 
distributions of these events are shown in Fig. (2). 



B. Effective signal lifetime 

We determine the effective lifetime of KgK^K~ decays from the distribution of 

proper decay times as follows. The proper decay time distribution of candidates can be 
parameterized as: 

Vit) =pIJ e~*'/" ■ Rit - t\ to)dt' + (1 - 'p)B{t), (10) 

where the first term represents the measured distribution of signal events with lifetime 
r, convolved with a resolution function, i?(t,to); corresponds to a possible shift of the 
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M„+„- [GeV/c2] Mk^k- [GeV/c2] 



FIG. 2: (Left) (M^+^-,M^o 

K+K-) distribution of selected events. (Right) Mj^+j^- distribution of 
events in the |M^o^+^- — m£)o| < 10 MeV/c^ and \M^+^- — mj^a \ < 10 MeV/c^ region (unfilled 
histogram), and 20 < |M^o^+^_ — m£)o\ < 30 MeV/c^ and iMj^+j^- — mj^o\ < 10 MeV/c^ region 
(filled histogram). Dashed (Bottom) vertical lines indicate the boundaries of ON (OFF) intervals. 



resolution function from zero, p = Ns/{Ns + A'^^) is the fraction of signal events, and the 
last term, B{t), describes the distribution of background events. Since the average of the 
convolution is the sum of the averages of the convolved functions, we can express the lifetime 
of signal events in region TZ (shifted for the resolution function offset) as 

Tn + to = , (11) 

where (t)'^ and (t)^ are the mean proper decay times of all events and background events, 
respectively. By measuring {t)^ and {t)^ for events in ON and OFF intervals of M^+k- we 
can obtain the two effective lifetimes and ycp from Eq. (9). Note that the resolution function 
offset, to, if small (to ^ t) and equal in ON and OFF regions, introduces a negligible bias 
(~ Vcp ' to/T) in the measurement, since it cancels in the numerator of Eq. (9). We use the 
simulated sample to confirm that the resolution function offsets t^^ and t^^^ are equal to 
within the statistical uncertainty. 

The requirement of minimal Kg candidate flight distance introduces a bias in the re- 
constructed mean proper decay time of signal KgK~^K~ decays: events where both 
D° and Kg candidates are short-lived are rejected by this requirement. This introduces an 
+0.5% bias in the mean of the measured proper decay times for D° KgK'^K'; the effect 
on the ycp parameter is smaller and is included in the systematic error. 



C. Signal and background fractions 

Signal and background fractions are determined from a fit to the distribution of events in 
the {Mj^+T^- , Mj^oj^+j^-) plane. In order to model the correlation between invariant masses 
^K^K+K- s-iid M,r+7r- of sigual eveuts (see Fig. 2), we parameterize the signal shape by a 
rotated triple two-dimensional Gaussian distribution. The individual Gaussians are required 
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FIG. 3: Invariant masses Mj^o^j^+p^^ (a, c) and M^+^- (b, d) of events passing all selection criteria 
for ON and OFF intervals in Mj^+j^-. Superimposed on the data (points with error bars) are 
results of the fit (solid line). 

to have the same mean value, which is allowed to vary in the fit. The ratio of the Gaussian 
widths is fixed to the Monte Carlo (MC) simulated value, and only the width of the core 
Gaussian and the three correlation coefficients are left free. 

Background events are classified into three categories according to their distribution in 
the (M^+TT-, M^o^+^-) plane (see Fig. 2): true Kg background, — > K^K'n'^ii" de- 
cays with the pion pair not originating from a Kg, and remaining background. True Kg 
background events are random combinations of charged kaons with correctly reconstructed 
Kg candidates; the shape in M^+^- is fixed to be the same as signal while in Mj^oj^+j^- it 
is parameterized with a second-degree polynomial. The remaining background events are 
random combinations of charged particles and are parameterized as a polynomial of first 
degree in M^+^- and second degree in M^o^+^-. The — >• K~^K~7i~^7i~ decays are peak- 
ing in Mj^oj^+j^-, but not in M^+^r-. According to MC simulation, the contribution of these 
events is small (~ 0.1%); thus they are not included in the fit but considered as a systematic 
uncertainty. 

The fractions and shapes are determined in a three-step fit for both ON and OFF regions. 
First, the fraction of signal events Fgig is obtained from a fit to the one-dimensional projection 
in Mj^oj^+j^-. In the second step, we fit the projection in M^+^- to find the sum of the 
fractions of signal and true Kg events, Fgig + -FtKs. Finally, we determine the signal shape 
parameters from a two-dimensional fit in which we use the Fgig and -FtKs results from the 
previous steps. The fitting procedure was checked using a high-statistics sample of simulated 



9 



25 
^ 20 
15 
10 
5 


-5 
-10 
-15 
-20 
-25 



; A 

: 


B 


- 






A 












1 
1 


IV 
1 V 










III 


- 








II 


























B 











-15 -10 -5 5 



10^ 



lOf^ 



10 



10 i; 



FIG. 4: The C) distribution of selected events. Signal box S, and sideband regions A and B are 
defined in the text. Quadrants denoted by I - IV are used in the systematic uncertainty estimate, 
as described in section IV. 



signal and background events and found to correctly reproduce the true event fractions. 

The results of this procedure are shown in Fig. 3. We find (72.3 ± 0.4) x 10^ signal 
events in the ON region and (62.3 ±0.7) x 10^ events in the OFF region. To achieve the best 
statistical accuracy on the ycp measurement, we optimize the size of the signal box. Because 
the invariant masses M^o^+^- and M^+^- are correlated for signal events, we define the 
signal box in the rotated variables: 

i = ^ (12) 

a/1 - V 1 ~ P ^° 

where Mko = 497.533 ± 0.005 MeV/c^ and M^o = 1864.874 ± 0.009 MeV/c^ are fitted K^g 
and masses, a^o = 1.880 ± 0.008 MeV/c^ and a^o = 2.839 ± 0.014 MeV/c^ are widths 
of the core Gaussian function, and p = 0.571 ± 0.003 is the correlation coefficient. The 
uncertainties are statistical only. The signal region that minimizes the statistical uncertainty 
on ycp (signal box) is found to be |^| < 3.9 and \(\ < 2.2. The two-dimensional distribution 
of (^, for the selected data is shown in Fig. 4. The signal fractions in the signal box are 
(96.94 ± 0.06)% and (90.53 ± 0.16)% in the ON and OFF intervals, respectively. 

The fraction of — > K~^K~7r~^7r~ decays in the signal box is estimated by fitting 
the Mj^oj^+j^- projection for events in the sideband regions M^+^- < 0.480 GeV/c^ and 
M^+^- > 0.514 GeV/c^, where the contributions of signal and true Kg background are 
small. The fractions of this background extrapolated to the signal box are found to be 
(0.02 ± 0.01)% and (0.07 ± 0.02)% in the ON and OFF intervals, respectively, and are 
reproduced well by MC simulation. 
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TABLE I: Mean proper decay times of events populating sideband regions A and B in the C) 
plane, {t)^ and (t)^, fractions (i = S, A, B) and estimated mean proper decay times of back- 
ground events, populating the signal box, for events in the ON and OFF intervals in Mj^+j^-. 
The uncertainties are statistical only. 



Mk+k- 


[fs] (t)B [fs] 


pS [%] pA [%] pB [%] 


{t)b N 


ON 
OFF 


223 ± 14 63.6 ± 4.7 
237.7 ±7.4 140.3 ±3.1 


93.31 ± 0.41 7.2 ± 1.8 91.83 ± 0.23 
90.17 ± 0.32 5.1 ± 1.2 88.02 ± 0.17 


60.8 ±4.8 
137.8 ±3.2 



D. Mean proper decay time of background events 

The mean proper decay time of background inside the signal box, {t)b, is determined from 
sideband regions A and B in the {C,, C) plane as shown in Fig. 4. The regions are chosen larger 
than the signal box to minimize the uncertainty on {t)h- To an excellent approximation, the 
mean proper decay times in sideband regions A and B ((t)^ and {t)^) can be expressed as 

(^)^ = PtKs(^)tKs±pit(^)rest, (14) 
(tf = PtKs(^)tKs±pLt(^)rest, (15) 

where and p^}^^ are the fractions of true Kg and the remaining background in region 
A (B). Similarly, the mean proper decay time of background in the signal box S can be 
expressed as 

{t)b=p!Ks{t)tKs+pLt{t)rcst (16) 

By solving Eqs. (14) and (15) for (t)tKs and (t)rcst5 and inserting the results into Eq. (16), 
we obtain 

\^)b = pA_pB ' y^') 

where = KKs/(PtKs + Plest)^ ^ = A,B,S. The fractions p*Ks and pj^^^, i = A, B, S are 
calculated from the results of the two-dimensional fit discussed in the previous section. In 
Table I we list the quantities used in Eq. (17) and the resulting {t)b for regions ON and 
OFF. 

In deriving Eq. (17), we have assumed that in regions A, B, and S the mean proper decay 
times (t)tKs and (t)rest are equal. This assumption has been validated using MC simulation. 
We have also neglected the signal leakage into regions A and B; if we compare, using MC 
simulation, the mean proper decay time of background events found in the signal box with 
that calculated from Eq. (17); we find agreement well within one standard deviation. The 
small deviations due to these assumptions are included in the systematic uncertainty. 



E. Fit to the sq distribution 

The Ai fractions, foN and foFF, are obtained from a fit to the Sq distribution. We use 
two different Dalitz models of K^K^K~ decays to parameterize the distribution: a 

four-resonance model from Ref. [21], and an eight-resonance model from Ref. [14]. The main 
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TABLE II: Fraction difference /on ~ /off for the two Dalitz models. The nominal values are 
calculated from the data in Refs. [14, 21], and the fitted values from our fit results. 





/on - 


/off 


Model 


Nominal 


Fitted 


4 res. [21] 


-0.730 ± 0.031 


-0.732 ± 0.002 


8 res. [14] 


-0.753 ± 0.004 


-0.769 ± 0.005 



sensitivity to ycp arises from J'^^0(1O2O) and i^^ao(980)° intermediate states, since the two 
have opposite CP eigenvalues. Because all resonance parameters cannot be determined from 
a one- dimensional fit, we fix the parameters of the resonances with smaller fit fractions using 
the amplitudes and phases from the corresponding model and world averages for masses and 
widths; we vary only the amplitudes of i^^0(lO2O) and i^'~ao(980)"'" (four-resonance model) 
or the amplitudes of i^50(lO2O) and i^~ao(1450)~'" (eight-resonance model), mass and width 
of the 0(1020), and the coupling constant qkk of the Flatte parameterization of the ao(980)°. 
The signal distribution is parameterized as 

^s(so)=e(so) j e{s+)\Ai{sQ,s+) + A2{sQ,s+)\^dsj,, (18) 

where e is the reconstruction efficiency determined from a sample of MC events in which the 
decay mode was generated according to phase space; the efficiency is found to be factorizable 
in the Dalitz variables sq and The background parameterization is obtained from the 
sideband region 5 < |C| < 25, where |^| < 3.9 corresponds to the signal region. A test of 
the MC So distributions of background events from the signal and sideband regions yields 
= 88.9 for 99 degrees of freedom; thus we conclude that the Sq distribution of events 
taken from the sideband region satisfactorily describes the background distribution in the 
signal box. 

Figure 5 shows fit results for the eight-resonance model, which we use to determine the 
fraction difference /on — /off, since it provides a better description of the Sq distribution. 
The reduced is 1-28 for the eight-resonance model and 1.91 for the four- resonance model 
for 230 degrees of freedom. In Table II the fraction differences /on — /off are given for both 
Dalitz models. The left column lists the values calculated from the data in Refs. [14, 21], 
and the right column lists the values calculated from the results of our fit. Uncertainties 
in /on — /off are calculated using the statistical errors of amplitudes and phases, without 
taking into account any correlation between them. Although the models are different, with 
distinct resonant structure [22], the differences /on ~ /off calculated for the two models are 
very similar. The small difference between them is included as a systematic uncertainty. 

F. Results 

Figure 6 shows the proper decay time distributions of selected events in the signal box 
S for the ON and OFF intervals. The distribution of background events is estimated from 
proper decay time distributions of events populating the sideband regions A and B and 
the known fractions of the true Kg background and the remaining background in all three 
regions. Inserting the values for (t), {t)b, and p (the fraction of signal) into Eq. (11) yields 
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1.0 1.05 1.1 1.15 1.2 



So [GeVVc^] 

FIG. 5: The sq distribution of — > KgK~^K~ decays with superimposed fit results for the 8- 
resonance Dalitz model given in Ref. [14]. The solid curve is the overall fitted function and the 
dashed curve represents the background contribution. 

TABLE III: Measured mean proper decay times in the signal box (t), effective background lifetimes 
(t)h, signal fractions p and the resulting effective signal lifetimes. The uncertainties are statistical 
only. 



Mk+k~ 


it) [fs] it), [fs] 


p[%] 


T + to [fs] 


ON 


402.7 ±2.5 60.8 ±4.8 


96.94 ±0.06 


413.4 ±2.5 


OFF 


386.7 ±2.6 137.8 ±3.2 


90.53 ±0.16 


412.7 ±3.0 



ToN+t^^ = (413.4±2.5) fs and ToFF+t^^^ = (412.7±3.0) fs. These results are summarized 
in Table III. The measured values for r ± to are close to the world average for t^o, and, 
since ycp <^ 1, this implies to/T is ~ 1% or less. Since the topology of events in the ON 
and OFF intervals is almost identical, we assume = t^^^ and include a systematic 
error to account for possible deviations from this assumption. This leads to a normalized 
lifetime difference (toff — ton)/(toff ± ton) = (—0.09 ± 0.47)% between the two regions. 
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t [ps] t [ps] 

FIG. 6: The proper decay time distributions of all events (unfilled histogram) and background 
events (hatched histogram) populating the signal box S for the ON and OFF intervals. 

where the uncertainty is statistical only. The difference in the Ai fraction corresponding to 
the eight-resonance model (see Table II) is foN — Ioff = (—0.769 ± 0.005); therefore, from 
Eq. (9) we obtain ycp = (+0.11 ± 0.61(stat.))%. 

IV. SYSTEMATICS 

We consider separately systematic uncertainties arising from experimental sources and 
from the D° — > K^K~^K~ decay model. First, we check the simulated sample to confirm 
that the resolution function offsets t^^ and t^^^ are equal. The small difference observed 
is consistent with the statistical error but conservatively propagated to ycp and taken as a 
systematic uncertainty (±0.38%). 

The mean proper decay time of background events populating the signal box (calculated 
from Eq. (17)) assumes a negligible contribution of signal events in sideband regions A and 
B, and also assumes equal mean proper decay times of the two background categories in all 
three regions. A, B and S. The systematic uncertainty resulting from the first assumption 
is evaluated by including the small residual fraction of signal events in regions A and B in 
the {t)b calculation; the resulting change in ycp is ±0.01%. The uncertainty resulting from 
the second assumption is evaluated by MC simulation; mean proper decay times of the two 
background categories are found to be consistent within statistical uncertainty in all three 
regions. Small differences between the mean proper decay times of the two background 
categories in the S, A and B regions result in ±0.09% and 0.04% variations of ycp for 
true Kg and remaining background, respectively. We add in quadrature the above three 
contributions to obtain a ±0.10% systematic error on ycp- 

The contribution of — > K~^K~7r^7r~ decays in our sample is found to be small and 
thus is not included. We evaluate their effect on ycp by taking the fraction of these events in 
the ON and OFF intervals from data, and their mean proper decay time from the simulated 
sample. The resulting change in ycp is ±0.07%. We include this change in the systematic 
uncertainty. 
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We study the choice of sideband regions used to determine {t)b as follows. The sidebands 
A and B are divided into four subregions (denoted I-IV) as shown in Fig. 4. The mean 
proper decay time of background events is then calculated using events in subregions (I, III) 
or (II, IV), and a difference of 0.05% in ycp is observed. This change is included as a 
systematic uncertainty. 

Possible systematic effects of selection criteria are studied by varying the signal box 
size and the selection criteria for at and the flight distance. Although no statistically 
significant deviation is observed, the maximum difference in ycp is (conservatively) assigned 
as a systematic uncertainty (±0.30%). 

The fitting procedure is tested using the simulated sample. A small difference between 
the fitted and true fractions of signal events in the signal box is propagated to ycp and 
included as a systematic uncertainty (±0.10%). 

The mean proper decay times of events populating the signal box S and the sideband 
regions A and B are taken to be the means of histograms of the proper decay times for events 
populating these regions. Changing the binning and intervals used in these histograms over a 
wide range results in a change in ycp of ±0.07%; we include this as an additional systematic 
uncertainty. 

Finally, we estimate the systematic uncertainty due to our choice of — > K^K^K~ de- 
cay model. First, we compare the fraction difference /on — /off obtained using the four- and 
eight-resonance Dalitz models. Despite the difference between the models in their resonant 
substructure [22], the values for /on — /off are similar (see Table II). We assign a 3% 
relative error to ycp due to the small difference in the above fractions. An additional 2% 
relative error is assigned due to the small difference between the fitted and nominal values of 
/on — /off- If the reconstruction efficiency e{s+) were constant, the contribution of the real 
and imaginary parts of the interference term A1A2 in Eq. (5) would vanish after integrating 
over s+. A slight decrease of e(s+) near the kinematic boundaries is observed from a large 
sample of simulated events; the effect of this variation on ycp is studied and found to be 
negligible. 

Adding all decay-model systematic uncertainties in quadrature with the statistical uncer- 
tainty in foN — foFF (= —0.769 ± 0.005, see Table II) yields a total uncertainty due to the 
decay model of 0.01%. Combining this in quadrature with all other sources of systematic 
uncertainty gives a total systematic error on ycp of 0.52%. The individual contributions to 
the total systematic error are listed in Table IV. 

V. SUMMARY 

We present the first measurement of ycp using a CP-odd final state in decays. Our 
method has the advantage of not requiring precise knowledge of the decay-time resolution 
function, and avoids several biases that can arise due to detector effects. The value of 
ycp obtained is 

ycp = (+0.11 ± 0.61(stat.) ± 0.52(syst.))%. 

This measurement of ycp using a CP-odd mode is consistent with previous measurements 
using CP-even final states [4, 5], and with the world average value ycp = (±1.13 ± 0.27)% 
[9]. 

We thank the KEKB group for the excellent operation of the accelerator, the KEK cryo- 
genics group for the efficient operation of the solenoid, and the KEK computer group and the 
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TABLE IV: Sources of the systematic uncertainty on ycp- 



Source 


Systematic error (%) 


Resolution function offset difference t^^^ — tg*^ 


±0.38 


Estimation of (t)^ 


±0.10 


K~7r~^7r" background 


±0.07 


Selection of sideband 


±0.05 


Variation of selection criteria 


±0.30 


Fitting procedure 


±0.10 


Proper decay time range and binning 


±0.07 


Dalitz model 


±0.01 


Total 


±0.52 
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APPENDIX: INTEGRATION OF AiA*2 OVER ONE DALITZ VARIABLE 

The amplitude A (A) for a (-0°) decay to a three-body final state, h^h~m^, depends 
on invariant masses of all possible pairs of final state particles: sq = M^+^_, s+ = M^+^o 
and S- = M^_^o- Only two of these three are independent, since energy and momentum 
conservation results in a constraint 

So ± s+ ± s_ = m|,o ± ml+ + ml- + m^o = m^. (A.l) 
In the limit of CP symmetry the following relation holds: 

A{so, s+) = A{so, S-) = A{so, w? - s+- So), (A. 2) 

and amplitudes Ai^2 (defined in section II) are then 

^i(so, s+) = ^ [^(so, s+) ± ^(sq, - s+ - So)] , (A. 3) 
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A2{so, s+) = ]- [A{so, s+) - A{so, - s+ - sq)] . (A.4) 



2 

We now show that 

/ Ai{so,s+)A*2{so,s+)ds+ = 0, (A. 5) 

where s™"(so) and s™'^^(so) are lower and upper bounds of Dahtz variable s+. For a given 
value of So, the range of s+ is determined by its values when the momentum of h'^ is parallel 
or antiparallel to the momentum of m°: 



v2 



where 



El. = (A.8) 

_ m^o - So + m^o 

- ^-7= lA.9j 



are the energies of and m° in the h^h rest frame. The left-hand side of Eq. (A. 5) yields 

-5™=''' (so) 

i(4-4 + J,- Jrf), (A.IO) 



/ A(so,s+)^;(so,s+)cis+ 



where 



/ ^(so,s+)^*(so,s+)ds+, (A.ll) 

^(so, s+)^*(so, — s+ — so)(is+, (A. 12) 



Ic= I A{so,w? - s+ - so)A*{sQ,s+)ds+, (A.13) 

Id= I A{sq, - s+ - sq)A*{sq, - s+ - so)ds+. (A.14) 

In integrals Ic and we perform a variable substitution s+ — > s_ (Eq. A.l): 

(is+ = —ds_, (A. 15a) 
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sT^'iso) sr(.o), (A.15b) 

sTiso) "'^^'^^ s^^%s,), (A.15c) 

and obtain 



Ic= A{sq, s^)A*{so,m^ - s_ - So)ds^ = (A.16) 

h= A{so,s^)A*{sQ,s.)ds. = la. (A.17) 



The right-hand side of Eq. (A. 10) therefore yields zero. 
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